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We reconcile the phenomenon of mesoscopic conductance fluctuations with the single parameter 
scaling theory of the Anderson transition. We calculate three averages of the conductance distri- 
bution: exp((ln<7)), (g) and 1/ (R) where g is the conductance in units of e 2 /h and R — 1/g is 
the resistance and demonstrate that these quantities obey single parameter scaling laws. We obtain 
consistent estimates of the critical exponent from the scaling of all these quantities. 



In the original proposal of the scaling theory of local- 
ization the conductance g (in units of e 2 /h) is the relevant 
parameter [jij. The scaling of this parameter is deduced 
by looking at the limiting cases of a good metal, where 
weak localization theory is applicable, and the strongly 
localized regime where an exponential dependence on sys- 
tem size is expected. A smooth monotonic interpolation 
is supposed between the two limits. The central equation 
of the theory is the j3 function 



ding 



(1) 



dlnL 

which describes how the conductance of an L x L x L 
cube is renormalized with system size L. 

The scaling theory has been very influential. It is the 
basis for the predictions that the Anderson transition is 
continuous and that the lower critical dimension for the 
Anderson transition is two. 

Shortly after its proposal, the pioneering numerical cal- 
culations of Pichard j| , and MacKinnon and Kramer || 
provided indirect support for the scaling theory. In these 
calculations the localization length of electrons on quasi- 
one dimensional bars was calculated. The theory of finite 
size scaling was then applied to deduce the critical pa- 
rameters from the dependence of the localization length 
on the transverse dimension of the bars. This has proved 
useful for making quantitative estimates of critical pa- 
rameters [@-||. 

It became clear later, in particular after the investiga- 
tion of universal conductance fluctuations in mesoscopic 
systems, that the conductance of a disordered system is 
a random variable |||l0| . In the vicinity of the Anderson 
transition the fluctuations in g are of the same order of 
magnitude as it's mean value ||. The smooth scaling 
behavior predicted by (Q) is clearly inconsistent with the 
fluctuating behavior of the conductance which occurs in 



practice and casts strong doubts over the soundness of 
the scaling theory ElJ]. There seem to be two principal 
remedies which we discuss in turn |l^,[l3| . 

First, we could attempt to establish that the distri- 
bution of conductance pl{9) (in the limit that the sys- 
tem size L and the correlation length £ are much longer 
than any microscopic lengths) obeys a single parameter 
scaling law. The precise meaning of this statement is 
that it should be possible to parameterize the bulk of the 
conductance distribution with a single parameter. If we 
denote this parameter by X then the bulk of the conduc- 
tance distribution must be of the form 



p L (g)~F(g;X(L)) 



(2) 



and the parameter X must obey a single parameter scal- 
ing law 



dlnX(L) 
dlnL 



f3(X) 



(3) 



Note that X need not necessarily be one of the moments 
of the distribution. Indeed, single parameter scaling of 
a distribution does not necessarily imply scaling of its 
moments. The moments may be dominated by non- uni- 
versal tails of the distribution or might not even exist. 
We shall refer to this first possibility as strong single pa- 
rameter scaling. 

Second, we could attempt to establish that some typ- 
ical or average conductance obeys a single parameter 
scaling law. This is a somewhat weaker statement since 
one parameter scaling of some average or typical quan- 
tity does not imply single parameter scaling of the cor- 
responding distribution. For example, it might happen 
that several independent parameters are needed to de- 
scribe the distribution. For this reason we shall refer to 
this as weak single parameter scaling. 
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Our purpose in this paper is to firmly establish single 
parameter scaling in the weaker sense given above, for 
systems near the critical point of the Anderson transition 
in three dimensions. A secondary objective is to identify 
which averages or typical values obey single parameter 
scaling laws. By simulating the conductance distribution 
for large ensembles of disordered Lx Lx L cubes we have 
achieved both objectives. 

We supposed perfect leads were attached to a pair of 
opposite sides of each cube and used the Landauer for- 
mula to relate the (Landauer) conductance <?l to the 
transmission matrix t which describes the transmission 
of electrons from one lead to the other llbl . 



g L = 2tr tt 



(4) 



We use the notation g^ to emphasize that this is the con- 
ductance that would be measured in a two probe mea- 
suring geometry. In the original work on the scaling the- 
ory the Thouless conductance g — Eq /A was considered 
where E c is the Thouless energy and A the mean energy 
level spacing. The Landauer conductance and the Thou- 
less conductance are not completely equivalent. Their 
relation has been considered in detail by Braun et al |l5| ] 
who state that the contact resistance, which is always 
present in a two terminal measurement, should be sub- 
tracted from gL ■ Hence we study the statistics of 



1 

f]L 



1 

2N 



(5) 



Here N = N(Ep) is the number of propagating chan- 
nels in the contacts at Fermi energy Ep and 1/2N is the 
contact resistance appropriate for the situation we have 
simulated. 

The motion of the electrons in the system is described 
by the Anderson model 



(6) 



<hj> 



where C|(Cj) denotes the creation (annihilation) oper- 
ator of an electron at the site j of a 3D cubic lattice. 
The amplitude of the random potential at site i is Wi. 
Hopping is restricted to nearest neighbors and its ampli- 
tude is taken as the unit of energy, V = 1. We supposed 
a box distribution with each Wi uniformly distributed 
on the interval [—W/2,W/2]. Previous work has veri- 
fied the universality of the critical behavior in this model 
with respect to the choice of distribution of the random 
potential P,E6 23 1. We imposed fixed boundary condi- 



tions in the transverse directions since we have found 
in previous work that corrections to scaling vanish more 
quickly with system size in this case j7|,|l7|. The Hamil- 
tonian (^|) has both time reversal and spin rotation sym- 
metries so that the observed critical behavior should be 
that of the orthogonal universality class. We used the 



method of Pendry et al |T|] to calculate the transmission 
matrix which appears in (3]). We set the Fermi energy 
Ep = Q.bV and for each combination of disorder W and 
system size L we generated an ensemble of 1,000,000 
samples (except for L — 16 where 500,000 were gener- 
ated). This allowed us to estimate the various averages 
of the conductance to roughly an accuracy of 0.1%. 

We examined the behavior of three different averages of 
the conductance distribution: exp ((In g)), (g) and 1/ (R) 
where R — l/g is the resistance. Note that for each av- 
erage it is possible to define a different (3 function. For 
example, for the typical conductance exp ((In <?}) 



/?(exp(0ns») 



d(ln.g) 
dlni ' 



(7) 



In the critical region the conductance fluctuations are of 
the same order of magnitude as the mean conductance 
so that these averages are not at all equivalent. 

We fitted the disorder and system size dependence of 
the averages to the standard scaling form. Taking the 
typical conductance as an example we supposed that 



(lug) = F (i)L 1/v ,(j)L y ^ 



(8) 



where tp is a relevant scaling variable and <p is an irrele- 
vant scaling variable which allows us to take account of 
corrections to scaling. Such deviations from perfect scal- 
ing are always present in a simulation of a finite system 
and it is necessary to have some means of accounting for 
them. We approximated this scaling function by its first 
order expansion in the irrelevant scaling variable and fit- 
ted the numerical data to the form 



(In 5 ) = F^L 1 '") + cj>V>F l ^L 1 ' l/ ). 



(9) 



The scaling variables were approximated by expansions 
in terms of the dimensionless disorder w = (W c — W) /W c 
where W c is the critical disorder separating the insulating 
and metallic phases. 



1p = IplW + 1p2W 



(10) 



The critical exponent v describes the divergence of the 
localization (correlation) near the critical point 



(ii) 



The absolute scale of the localization length £o cannot 
be determined from this fit. The system size dependence 
of the irrelevant scaling variable is described by an ex- 
ponent y < 0. The functions Fq and F± were expanded 
to third order in w. We found that this fitting scheme 
was the simplest which still allowed for goodness of fit 
probabilities in excess of 0.1. 

The results of the analysis are displayed in Table | and 
Figures [M. A number of points can be noticed. First, 
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acceptable fits to the single parameter scaling law are ob- 
tained for all three averages. Second, mutually consistent 
values for the critical exponent are obtained and these es- 
timates are also consistent with estimates based on the 
scaling of the localization length of electrons on bars || . 
Third, mutually consistent values of the critical disorder 
are obtained. Finally, the critical values of the quantities 
considered vary widely, indicative of the broadness of the 
distribution of the conductance fluctuations in the crit- 
ical regime. This means that the precise form of the j3 
function is different for each average, though all should 
have the same slope (= \jv) at their respective zeros 
|Q| . Note that the /3 function will also depend on the 
boundary conditions since it is known that critical con- 
ductance distribution depends strongly on the boundary 
conditions even in the limit that Z/ — > oo [Q] . 

These results firmly establish that the typical conduc- 
tance, the mean conductance and the mean resistance all 
obey single parameter scaling in the critical regime. 

To what extent do these results also support scaling 
in the strong sense? Let us consider the vicinity of the 
metallic, insulating and critical fixed points in turn. 

The metallic fixed point can be reached by taking the 
limit L — > oo with £ fixed from any starting point on the 
metallic side of the transition. In this limit PL{g) ap- 
proaches a normal distribution with a size independent 
variance 12 1. Therefore a single parameter, (g), is suf- 
ficient to parameterize the distribution. Further it can 
be established that (g) obeys a single parameter scaling 
law using weak localization theory near the metallic fixed 
point @. 

The insulating fixed point can be reached by taking the 
limit L — » oo with £ fixed from any starting point on the 
insulating side of the transition. In the insulating regime 
exact results are available only for one dimensional sys- 
tems [^9|. It is found that distribution of conductance is 
log- normal with a width which is related to the mean 
of the distribution (provided certain conditions are satis- 
fied pp[.) Hence a single parameter (In (7) is sufficient to 
parameterize the distribution. Also (In g) obeys a single 
parameter scaling law in the limit L ^> £. Calculations 
of the average conductance and the average resistance, 
in this limit, in strictly one dimension, give 



of the occurrence of rare resonance states with anoma- 
lously large conductances which are not described by the 
the log- normal distribution |l9|. It is possible, but not 
certain, that these results are qualitatively correct for the 
insulating fixed point in three dimensions. 

At the critical fixed point £ diverges and the sin- 
gle parameter scaling hypothesis predicts a scale in- 
dependent universal critical conductance distribution 
p c {g)- This has been confirmed in numerical simulations 

Thus, it seems likely that the conductance distribution 
also obeys single parameter scaling in the strong, as well 
as the weak, senses. However, the available results are 
not completely conclusive. It seems to us that a demon- 
stration that Ph{g) obeys (^|) and (§) in the limit L — > 00 
and £ — > 00 with L/£ fixed for any value of the ratio L/£ 
is required. 
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X 


exp((lng)) 


(9) 


l/(R) 


V 


1.57(56,58) 


1.58(57,60) 


1.54(53,56) 


w c 


16.48(47,49) 


16.47(45,48) 


16.49(48,50) 


x c 


0.291(290,293) 


0.573(570,576) 


0.100(099,101) 


Q 


0.5 


0.4 


0.3 



TABLE I. Results of the scaling analysis for each average 
X. Estimates of the critical exponent v and the critical disor- 
der Wc together with 95% confidence intervals are given. Xc 
is the value of the relevant statistic at the critical point. The 
number of data values was 199 and the number of parameters 
in each fit was 12. The goodness of fit probability Q is also 
given. 



(fj) 



2L 



3/2 



cxp (-L/2H) 



(12) 



and 



(R) = - exp (4L/0 



(13) 



showing that both these quantities also obey single pa- 
rameter scaling laws. However, the average conductance 
(and also its higher moments) are not consistent with the 
log- normal form of the bulk of the conductance distribu- 
tion. They turn out to be determined by the probability 
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FIG. 1. The logarithm of the typical conductance versus 
the amplitude of the potential fluctuations for system sizes 
L = 4, 6, 8, 10, 12, 14 and 16. The solid lines are the fit to the 
data. 
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FIG. 2. The same data as in Fig. 1 after corrections to 
scaling are subtracted and plotted versus L/£ to exhibit the 
single parameter scaling function. 
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